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Abstract 



. Kaleu is an independent, true phase space generator. After providing it with 

some information about the field theory and the particular multi-particle scatter- 
ing process under consideration, it returns importance sampled random phase space 
■^j- ■ points. Providing it also with the total weight of each generated phase space point, it 

^ \ further adapts to the integration problem on the fly. It is written in FORTRAN, such 

O that it can independently deal with several scattering processes in parallel. 

o ' 



1 Introduction 



Efficient phase space generation is an important issue in the study of multi-particle processes at 
collider experiments. Parton-level defined observables are eventually formulated as phase space 
integrals, and the evaluation of these is the most time-consuming part of a calculation. The 
complicated structure of both the integrand and massive multi-particle phase space enforces the 
use of Monte-Carlo methods, setting the total evaluation time to the sum of the generation of 
phase space points and the evaluation of the integrand at the phase space points. While it is 
important to keep the computing time of each integrand evaluation as low as possible, it is also 
important to keep the number of evaluations to reach an acceptable accuracy as low as possible. 
The latter is commonly referred to as efficient phase space generation, and eventually is a trade 
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off between the number of evaluations and the average computing time for a single generation of 
a phase space point. 

Since efficient phase space generation should already be settled at leading order (LO) cal- 
culations, automatic systems for LO calculations developed over the last decade deal with this 
issue [|2]-riTi. Various techniques are used in various combinations. Efficient phase space gen- 
eration implies knowledge about the integrand, which is utilized both before the integral is per- 
formed, through phase space mappings, and during the integration process, using the actual value 
of the integrand at the phase space points. Methods using the latter are referred to as adaptive 
methods, and the ones used within particle physics are based on adaptive grids dividing phase 
space into subspaces 1(72141611 . concern the optimization of the parameters in mixture distribu- 
tions (the so-called adaptive multi-channel method) IfTTll . or combine both lfT8ll . Most important 
is, however, to use as much as possible information about the integrand before performing the 
integral, be it just when deciding which underlying phase space mapping to use in adaptive grid 
methods. 

Eventually, the aim of the game is to have the probability density with which phase space 
points are generated look as much as possible as the integrand and to let it have the same peak 
structure. The fact that, for parton-level multi-particle calculations, this structure is essentially 
determined by the sum of the squares of the tree-level Feynman graphs contributing to the scat- 
tering amplitude is of great help. Since each Feynman graph encodes the phase space mapping 
needed for the efficient integration of its own square, the multi-channel method, in which each 
graph represents a channel, can be applied for efficient performance of the full integration. The 
problem with this approach is that the number of Feynman graphs grows factorially with the 
number of final-state particles. This problem was solved for the calculation of the tree-level 
amplitudes themselves through the use of recursive techniques [|8l [T9] - l23l . Regarding the phase 
space generation, the straightforward solution within the multi-channel method is to keep only a 
restricted number of channels with the highest channel weight. However, in order to determine 
these, the integration process still has to be started including all channels. 

The problem only exists in the calculation of the weight coming with each phase space point. 
Only there the sum over all graphs occurs, not in the actual generation of the phase space points, 
for each of which only one graph encoding a phase space mapping is chosen. The reason why the 
weight cannot be calculated with the same recursive approach like the amplitude is that, within 
the multi-channel method, each individual graph has its own unique weight. And indeed, the 
problem is solved by departing from this approach, as was pointed out in [|25l . and instead giving 
individual weights to the vertices. Instead of choosing a whole graph encoding a mapping at 
once for each generation of a phase space point, a mapping is composed by choosing respective 
branchings. Each branching corresponds to a vertex, and has its unique weight among the list of 
possible branchings at a given stage in the construction of the phase space mapping. The weight 
associated with the generated phase space points can then be calculated following recursion 
relations analogous to those for the amplitude calculation. 

As said, the foregoing works well assuming the peak structure of the integrand is determined 
by sum of the squares of the tree-level Feynman graphs contributing to the scattering amplitude, 
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which however is not a gauge-invariant quantity. For processes involving many gluons, other 
phase space mappings have been designed ||26] - |28l . 

Although many implementations of the various methods for efficient phase space generation 
exist, most of them are heavily integrated in automatic systems for full multi-particle calcula- 
tions, making it difficult to use them separately. This write-up presents a program with which 
Kinematics Are Lucidly and Efficiently Utilized^. It makes a number of routines available which 
can be used for phase space generation in any given Monte Carlo program. As part of the ini- 
tialization, the user has to provide masses and widths of all possible real and virtual particles in 
the multi-particle process under consideration, and a list of interaction vertices. Kaleu uses this 
information to efficiently map phase space following the method of the Recursive Phase space 
Generator described in Il25ll . The generator is optimized on the fly during the integration pro- 
cess, with the option to also optimize the generation of invariants and polar angles further with 
Parni |[T6ll . It is a Fortran program, written such that several instances of the program can 
operate in parallel, completely independently of each other. It can be obtained from 

|http : //annapurna . if j . edu . pi/ "hameren 

2 The algorithms 

As mentioned before, KALEU uses the algorithm of the Recursive Phase space Generator from 
11251 . Like most algorithms, it uses the fact that n-body phase space can be decomposed in a 
sequence of two-body phase spaces, referred to as splittings in the following. 

2.1 Two-body phase space generation 

In each splitting, two new four-momenta pi,p2 are generated that sum up to an, possibly pre- 
viously generated, existing four-momentum Q = p! + p 2 . Depending on whether 0, 1 or 2 of 
these momenta are on-shell, 2, 3 or 4 random variables have to be generated from which p! , p 2 
are constructed. Most common is to generate the invariants pf, p§, the cosine z of the polar angle 
of pi in the center-of-mass frame (CMF) of Q, and the the azimuthal angel c() in the same frame. 
If any of the momenta is on-shell, and less than 4 random variables have to be generated, we say 
"the invariant is generated with a Dirac delta-distribution". 

The azimuthal angle is usually not correlated with any existing kinematical variable, but the z- 
variable sometimes is. It may be the cosine of the polar angle w.r.t. to another, already generated, 
momentum in the CMF of Q. The z- variable is then linearly related to the scalar product of 
the two momenta. This happens, for example, in the so called open antenna generation in [|28l 
in order to arrive at antenna densities for the efficient integration of multi-gluon amplitudes. It 
also happens in the the so-called t-type generation UJ. In this case, the z- variable is taken to 
be the cosine of the polar angle w.r.t. to one of the initial-state momenta qi in the CMF of Q, 
and is linearly related to t = (pi — q! ) 2 = (Q — qi — P2) 2 . If this generation is formulated as 

1 Kaleu outranks Sarge. 
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Table 1 : List of vertices 
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for the process uu — > ddZ 



(qi, Q') — > (pi,p2) with Q' = Q — qi, these four momenta can be interpreted as the external 
legs of a t-channel Feynman graph. This is useful for bookkeeping purposes, as it will also be 
for Kaleu. Both this t-type generation and the "s-type" generation with uncorrelated z variable 
are applied in KALEU. 

The mentioned different choices of random variables are applied in order to have control 
over the densities following which they are distributed. They should match the behavior of the 
the squared matrix element one is trying to integrate. KALEU uses the densities as mentioned for 
example in H|. 



2.2 Recursive phase space decomposition 

We will illustrate how sequences of splittings are composed with the help of an explicit ex- 
ample, namely the process uu — > ddZ. For this process, Kaleu will create the list of split- 
tings/mergings, or just vertices, in Tabled] The list contains almost all possible three-point ver- 
tices occurring in Feynman graphs for this process, i.e. all three-point vertices to be calculated 
in a recursive calculation of the amplitude. A few are missing, as will be explained below, and a 
few appear several times with differences in the last column, which will also be explained. The 
list only contains three-point vertices, the meaning of the particles between square brackets is 
explained below. There is no distinction between particles that only differ in charge since it does 
not matter for the kinematics. The number between parentheses encodes the momentum in the 
binary representation. The external particles have momenta u(l )u(l 6) — > d(2)d(4)Z(8), and 
for example the gluon/photon/Z-boson attached to the d(2)d(4) pair has momentum 2 + 4 = 6. 
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For a process with n particles in the final state, momentum conservation dictates that initial-state 
momentum 2 n+1 is equal to minus momentum 2 n+1 — 1 =1 + 2 + 4 + . .. + 2 n . 

To generate a phase space point, the list should be read upside down, and each vertex rep- 
resents a splitting. The starting point is the sum of all final-state momenta 2 + 4 + 8, plus one 
initial-state momentum 1 . The latter is necessary to include t-type kinematics. Thus a vertex is 
chosen among number 18 to 33. Each of these carries its own multi-channel weight, which is 
updated during the Monte Carlo process. Suppose vertex 29 is chosen, having W(3) , d(12) on 
the r.h.s. Then the same is repeated for all vertices with W(3) on the l.h.s. and with d(12) on 
the l.h.s.. For the latter, this involves only vertex 5. For the former, no such vertex exists in the 
list, because it would represent the trivial operation of subtracting the initial-state momentum 1 
from 3 to obtain the final-state momentum 2. Thus the constructed Feynman graph representing 
the phase space mapping is given by 



u(15) 




(1) 



Vertex 29 represents a t-type splitting, as all vertices do with an odd momentum on the l.h.s., 
because they contain momentum 1 . The first step in a t-type splitting is the generation of the 
positive invariants, in the example above the squares of momenta 12 and 2. The latter is actually 
external and fixed to the squared mass of the d-quark, so it is "generated following a delta- 
distribution". One of these invariants, in this case the square of momentum 2, is beyond the 
encoding of the three-point vertex, and has to be listed separately. This is the purpose of the 
particles between square brackets in the list of vertices. Splittings of the s-type, like vertex 12, 
do not have entries in this last column, as well as t-type vertices which give the possibility to be 
followed by another t-type splitting, like vertex 25. 



u(15) 



(2) 



W(3) \ 
u(1) d(2) 



As mentioned before, vertices representing a splitting involving momentum 1 explicitly are omit- 
ted. The only exception are vertices with the sum of all final-state momenta on the other leg, i.e. 
vertices 1 8, 1 9 and 20, since they are the starting point for pure s-type Feynman graphs like 



u(15) 




(3) 
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For the calculation of the weight associated with a phase space pointed generated as above, 
the list of vertices has to be read in the order as it stands. Each "< — >" should now be inter- 
preted as "= +" i.e. each vertex represents a contribution to the "off-shell current" labeled by 
particle-type and momentum on the l.h.s.. These off-shell currents are probability densities, and 
the final density labeled by u(15), obtained after executing the whole list of vertices, is the re- 
ciprocal of the weight. Each vertex includes the density associated with the variables generated 
for the splitting, and the multi-channel weight for the choice of the vertex among the possibilities 
contributing to the same off-shell current. This way, the contribution of all Feynman graphs is 
included, underlining the recursive character of the algorithm. 

3 Usage 

Kaleu has been designed such that the user can conveniently write his own user-defined routines 
suitable for a given Monte Carlo program. Some examples are included in the program. In 
particular, an interface to replace Phegas with Kaleu in the Helac/Phegas system RH23II24H 
is included. 

The first thing the user has to do is provide Kaleu with all particles and possible vertices. 
For the example in Section |2] the following routine would be sufficient: 

subroutine my_model ( model ,vertx ) 
use avh_kaleu_model 

type (model_type ) , intent (out) model 
type ( vertx_type ) , intent (out) :: vertx 
integer :: gluon, photon, wboson, zboson, uquark, dquark 

j 

parameter) gluon=l ) 
parameter ( photon=2 ) 
parameter ( wboson=3 ) 
parameter ( zboson=4 ) 
parameter) uquark=5 ) 
parameter ( dquark=6 ) 



call 


addparticle ( 


model 


, gluon 


,'g ' 


, OdO 


, OdO 




call 


addparticle ( 


model 


, photon 


, ' A ' 


, OdO 


, OdO 




call 


addparticle ( 


model 


, wboson 


, 'W ' 


, 80 . 419d0 


,2.04 


3d0 ) 


call 


addparticle ( 


model 


, zboson 


, ' Z ' 


, 91 . 188d0 


,2.44 


5d0 ) 


call 


addparticle ( 


model 


, uquark 


r ' U ' 


, OdO 


, OdO 




call 


addparticle ( 


model 


, dquark 


, ' d ' 


, OdO 


, OdO 





call add/vertex ( vertx , uquark, uquark, gluon ) 

call addvertex ( vertx , dquark, dquark, gluon ) 

call addvertex ( vertx , uquark, uquark, photon 

call addvertex ( vertx , dquark, dquark, photon 

call addvertex ( vertx , uquark, uquark, zboson 

call addvertex ( vertx , dquark, dquark, zboson 

call addvertex ( vertx , uquark, dquark, wboson 

call addvertex ( vertx , wboson, wboson, zboson 

call addvertex ( vertx , wboson, wboson, photon 

call addvertex ( vertx , gluon, gluon, gluon ) 

end subroutine 
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The output model and vertx are of public derived type with private components, so the user 
can only declare them and pass them to routines from Kaleu, in this case addpart icle and 
addvertex. The input of these routines should be clear, it should only be mentioned that the 
order of the particles in addvertex does not matter. In practise, the user would of course 
extend the model to, for example, the full standard model, which was refrained from for brevity 
here. Some prepared model routines are included in the program. 

A routine as above would typically be called in a user-defined initialization routine to be 
called in the Monte Carlo program. A complete set of such user-defined routines could look as 
follows: 

module my_kaleu 

use avh_kaleu 

use avh_kaleu_model 

use avh_kaleu_kinem 

type (model_type) , save 

type (kinem_type) , save 

type (kaleu_type) , save 
end module 



mdl ! Masses and widths 
kin ! Kinematics 

obj ! Instance of the phase space generator 



subroutine my_init ( process , nfinst , ecm , nbatch, nstep, thrs ) 
use my_kaleu 

integer , intent (in) :: process ( -2 : 1 7 ), nfinst , nbatch, nstep 

real (kind ( IdO ) ) , intent (in) :: ecm, thrs 
type (vert x_type) :: vtx 
call my_model ( mdl, vtx ) 

call kaleu_put_process ( mdl, vtx, kin, obj , process , nfinst ,ecm ) 
call kaleu_init_adapt ( mdl, obj , nbatch, nstep, thrs ) 
end subroutine 

subroutine my_gnrt ( discard , pkaleu ) 
use my_kaleu 

logical , intent (out) :: discard 

real (kind (IdO) ) , intent (out) :: pkaleu ( : 3 , -2 : 1 7 ) 
call kaleu_gnrt ( mdl, kin, obj , discard , pkaleu ) 
end subroutine 

subroutine my_wght ( weight ) 
use my_kaleu 

real (kind ( IdO ) ) , intent (out) :: weight 
call kaleu_wght ( mdl, kin, obj , weight ) 
end subroutine 

subroutine my_collect ( weight ) 
use my_kaleu 

real (kind ( IdO ) ) , intent (in) :: weight 
call kaleu_collect ( obj , weight ) 
end subroutine 

The module my_kaleu contains essentially what in a traditional Fortran program would 
be the common blocks. The items are of public derived type with private components, so the 
user can only declare them and pass them to routines from Kaleu. By declaring them as arrays, 
several instances of the program can be created. The user-defined routines above are wrappers 
of Kaleu -routines which need the items in my_module as arguments. 

The initialisation routine my_init takes the process as input, and the particles should 
have the same integer labelling as in the user-defined my_model. The entries —2 and —1 of 
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process are the initial-state particles, the (strictly) positive entries contain the final-state par- 
ticles. The integer nf inst is the number of final- state particles. The input ecm is the center- 
of-mass energy. The numbers nbatch, nstep and thrs are related to the multi-channel 
optimization. This happens in nstep steps of nbatch non-zero weight events. After these 
nstep steps, each channel with a weight smaller than thrs times the average channel weight 
at the off-shell current it belongs to is discarded. 

Kaleu determines limits on two-particle invariants based on the masses of final- state par- 
ticles, but these may be equal zero and in a practical calculation there will be phase space cuts 
present. The user may translate these into two-particle invariant mass cuts, and feed them to 
Kaleu for more efficient phase space generation with 

call kinem_updt_smin ( kin , smin ) 

The entries smin ( 1 : nf inst , 1 : nf inst ) of the array smin (-2:17,-2:17) should be 
minima to final-final state two-particle invariants, and the entries smin (-2 : -1 , 1 : nf inst ) 
should be negative, and should be maxima to final-initial state two-particle invariants. The array 
is considered to be symmetric. 

All arguments of the generation routine my_gnrt are output. If the logical discard is 
true, the phase space point should get weight zero, and the array pkaleu must not be used. If 
discard is false, pkaleu contains the momenta, including those of the initial-state particles. 
The third routine my_wght returns the weight associated with the most recently generated phase 
space point, and the final routine my_collect takes the full weight, including the integrand, 
as input for optimization purposes. 

The routines above would suffice to deal with e + e scattering. In case the user would like to 
generate events for hadron-hadron scattering, initial-state momenta can be provided to KALEU 
for each event with 

call kinem_inst ( kin , discard , pkaleu ) 

This routine should be called just before kaleu_gnrt. The array pkaleu has the same for- 
mat as before, but is input here. The entries pkaleu (0:3,-2) and pkaleu (0:3,-1) are 
taken as the initial-state momenta. The logical discard is output, and is true if something 
is wrong with the momenta. Alternatively, the user can let Kaleu deal with the x-variables 
within collinear factorization. For this option, one item should be added to the module, and a 
few routine calls should be added. An example is given in Appendix [A] The optimization of the 
generation of the x-variables is performed with the help of Parni. 

An example of the use of more instances of the program is given in Appendix [Bj 

4 Results 

In order to assess the performance and applicability of KALEU, some results existing in literature 
are reproduced. 
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Table 2: Cross sections in [fb] for some e + e —> 6f processes as in [29]. The numbers in the 
column "AMEGIC++" and "H/PHEGAS" are taken directly from that paper. The number N eval 
refers to the run with "H/Kaleu". 

Tabled Tabled Table|4]and Table[5]show some results for e + e~ — > 6f processes as presented 
in [|29l . It concerns all processes at y/s = 500 GeV, and "with QCD" if applicable. The numbers 
in the columns "AMEGIC++" and "H/Phegas" are taken from that paper. The numbers in the 
column "H/Kaleu" were obtained by replacing Phegas with Kaleu in the Helac/Phegas 
system. In [|29l results where obtained using 10 6 phase space points before cuts. It is, however, 
not clear whether this includes an optimization phase. The results with KALEU were obtained 
with 1 7 phase space points before cuts, including optimization. In order to compare the perfor- 
mances, both the standard deviation o and vTO cr are presented explicitly for the latter. Some 
results with "H/Phegas" have been rounded further than in the original paper for clarity. This 
leads to an ambiguity for one value, which has been highlighted. The optimization was per- 
formed in 10 steps of 50 x 10 3 phase space points after cuts. These were not included in the 
estimation of the cross section. The last column in the tables presents the total number of phase 
space points passing the cuts and leading to a matrix element evaluation, including the 500 x 1 3 
evaluations used for optimization. 

Table |6] shows some results for e + e~ — > 8f processes as presented in QUI . It concerns 
some processes with iuh = 130 GeV and m t = 1 74.3 GeV, at y/s = 800 GeV and with QCD. 
The numbers in the column "Carlomat" are directly taken from that paper. The numbers in 
the column "Helac/Kaleu" were again obtained by replacing Phegas with KALEU in the 
Helac/Phegas system. They were obtained with 1 8 phase space points before cuts, including 
optimization. The optimization was performed in 1 steps of 1 00 x 1 3 phase space points after 
cuts, which were not included in the estimation of the cross section. The last column in the tables 
presents the total number of phase space points passing the cuts and leading to a matrix element 
evaluation. 

In order to present some results for hadron-hadron collisions, Kaleu has been connected 
with Alpgen [7J. This program deals with partonic subprocesses in hadron-hadron collisions 
in one Monte Carlo run, making an optimized choice of a subprocess to be generated for each 
event. Because Kaleu can provide an independent generator for each subprocess, it can easily 
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Higgs 


■v e v e ud du 


0.4755(21) 


0.4711(24) 


0.4745(1 8)(06) 


7.135e+6 




y e y e ude~-v e 


0.16033(63) 


0.16011(78) 


0.16123(47)(15) 


9.018e+6 




y e y e ud u. _ v H 


0.14383(53) 


0.14439(65) 


0.14367(36)(11) 


9.239e+6 



Table 3: Cross sections in [fb] for some e + e — > 6f processes as in [29J. The numbers in the 
column "AMEGIC++" and "H/Phegas" are taken directly from that paper. The number N eva i 
refers to the run with "H/Kaleu". 





Final state 


AMEGIC++ 


H/Phegas 


H/Kaleu(\/10o-)(1o-) 


Neval 










0.03747(29) 


0.03749(32) 


0.03740(3 1)( 10) 


4.868e+6 


with 






ud e~v e 


0.1106(22) 


0.1090(07) 


0.1083(10)(03) 


5.245e+6 


Higgs 






e~e + 


2.731(065)e-3 


2.691(042)e-3 


2.737(108)(034)e-3 


2.285e+6 








uudd 


0.2634(22) 


0.2642(15) 


0.2638(22)(07) 


5.218e+6 








uuuu 


8.767(65)e-3 


8.978(58)e-3 


8.882(87)(28)e-3 


3.614e+6 






> 4 




0.03054(23) 


0.03092(19) 


0.03079(22)(07) 


4.512e+6 


no 




> 4 


ud e~v e 


0.08911(53) 


0.08925(48) 


0.08932(51)(16) 


4.800e+6 


Higgs 






u~u. + e~e + 


2.280(66)e-3 


2.277(62)e-3 


2.224(82)(26)e-3 


1.847e+6 








uudd 


0.2092(12) 


0.2075(13) 


0.2092(1 0X03) 


5.152e+6 






> 4 


uuuu 


6.134(29)e-3 


6.108(27)e-3 


6.072(41)(13)e-3 


3.907e+6 



Table 4: Cross sections in [fb] for some e + e — ) 6f processes as in [29 1. The numbers in the 
column "AMEGIC++" and "H/Phegas" are taken directly from that paper. The number N eval 
refers to the run with "H/Kaleu". 
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Final state 



AMEGIC++ H/Phegas H/Kaleu(vT0o-)(1ct) 



N 



eval 



withHiggs fi-ji+bbbb 30.96(0.60)e-3 30.19(0.43)e-3 29.96(1. 42)(0.45)e-3 3.153e+6 



noHiggs u~u + bbbb 6.308(24)e-3 6.364(21)e-3 



6.377(26)(08)e-3 4.918e+6 



Table 5: Cross sections in [fb] for some e + e 6f processes as in [29]. The numbers in the 
column "AMEGIC++" and "H/Phegas" are taken directly from that paper. The number N eval 
refers to the run with "H/Kaleu". 



Final state 


Carlomat 


Helac/Kaleu 


^eval 


bb bb v^u. + t~v t 


35.9(1) 


36.17(15) 


18.12e+6 


bb bb v e e + e~y e 


36.1(1) 


36.37(18) 


16.82e+6 


bb bb ud yTv^ 


100.6(2) 


100.50(46) 


13.21e+6 


bbbb cs e~v e 


100.5(3) 


100.52(47) 


12.63e+6 


bb bb ud sc 


314(2) 


314.1(1.4) 


13.50e+6 


bb bb ud du 


314(1) 


320.1(1.9) 


13.36e+6 



Table 6: Cross sections in [ab] for some e + e — > 8f processes as in [|30l . The numbers in the 
column "Carlomat" are taken directly from that paper. The number N eva i refers to the run with 
"Helac/Kaleu". 

be merged into the structure of Alpgen. This has been done such, that Alpgen still makes 
the choice of a subprocess for each event, and that, given the subprocess, Kaleu generates the 
variables Xi , Xi for the PDFs and the final- state momenta. It should be noted, however, that for 
the purpose of this write-up, the connection between the two programs has been established only 
to the level of cross section calculation and not, for example, to full event un- weighting. Table [7] 
presents results for processes ofthetypepp — > tt+Njets and pp — > ttbb + Njets. The numbers 
in the column "Alpgen" are taken directly from [7J. The user of Alpgen chooses the number 
of matrix element evaluations for a Monte Carlo run, not the number of generated phase space 
points before cuts. The last column presents the number of evaluations for the run with Kaleu, 
including those used for optimization but omitted for the estimation of the cross section. 



5 Summary 

The program Kaleu for parton-level phase space generation was presented. Given the masses 
and widths of all real and virtual particles, it generates importance sampled phase space points 
for the scattering process under consideration. Given the total weight of each event, including the 
value of the integrand, it optimizes further to this integrand on the fly. It is written in Fortran, 
such that it can independently deal with several scattering processes in parallel. 



11 



Final state 


Alpgen 


Alpgen/Kaleu 


Neval 


tt + 2jets 


255(1) 


254.38(73) 


11 A „ i ZT 

1 1.4e+6 


tt + 3jets 


111.5(5) 


1 1 1 C\C\ S AO\ 

111.09(43) 


11.4e+6 


tt + 4jets 


A O A / A\ 


41. /2(45) 


1 1 £~ i ^ 
1 1.6e+6 


tt + ojets 


1 A f\H / I £ \ 

14.U/(16) 


14.36(13) 


1 1 £~ i ^ 
1 1.6e+6 


tt + 6jets 


4.36(8) 


4.369(43) 


13.0e+6 


ttbb 


1.35(1) 


1.3490(21) 


10.2e+6 


ttbb + Ijet 


1.47(2) 


1.4624(42) 


10.4e+6 


tt bb + 2jets 


0.94(2) 


0.9280(40) 


10.6e+6 


tt bb + 3jets 


0.457(8) 


0.4522(28) 


10.6e+6 


ttbb + 4jets 


0.189(4) 


0.1851(14) 


6.2e+6 



Table 7: Cross sections in [pb] for some pp collision processes at LHC as in [J7| . The numbers 
in the column "Alpgen" are taken directly from that paper. The number N eval refers to the run 
with "Alpgen/Kaleu". 
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A Generation of x- variables 



An example of user-defined routines which can be used for the generation of events for hadron- 
hadron scattering within collinear factorization. 

module my_kaleu 

use avh_kaleu 

use avh_kaleu_model 

use avh_kaleu_kinem 

type (model_type ) , save :: mdl 

type (kinem_type) , save :: kin 

type (kaleu_type) , save :: obj 

type (strf_type) , save :: str 
end module 

subroutine my_init ( process , nf inst , ecm , nbatch, nstep, thrs ) 
use my_kaleu 

integer ,intent(in) :: process (-2 : 17 ), nf inst , nbatch, nstep 

real (kind ( IdO ) ) , intent (in) :: ecm, thrs 
type (vert x_type ) :: vtx 
call my_model ( mdl, vtx ) 

call kaleu_put_process ( mdl, vtx, kin, obj , process , nfinst , ecm ) 
call kaleu_init_adapt ( mdl, obj , nbatch, nstep, thrs ) 
call kaleu_init_strf ( str, kin , OdO ) 
end subroutine 

subroutine my_gnrt ( discard , xlkaleu, x2kaleu, pkaleu ) 
use my_kaleu 

logical , intent (out) :: discard 

real (kind ( IdO ) ) , intent (out) :: , xlkaleu, x2kaleu, pkaleu ( : 3 , -2 : 17 ) 
call kaleu_gnrt_strf ( str, kin , discard , xlkaleu, x2kaleu ) 
call kaleu_gnrt ( mdl, kin, obj , discard , pkaleu ) 
end subroutine 

subroutine my_wght ( weight ) 
use my_kaleu 

real (kind ( IdO ) ) , intent (out) :: weight 
real (kind(ldO) ) : : ww 

call kaleu_wght ( mdl, kin, obj , weight ) 
call kaleu_wght_strf ( str , ww ) 
weight = weight*ww 
end subroutine 

subroutine my_collect ( weight ) 
use my_kaleu 

real (kind ( IdO ) ) , intent (in) :: weight 
call kaleu_collect ( obj , weight ) 
call kaleu_collect_strf ( str , weight ) 
end subroutine 

The last argument of kaleu_init_st rf is a minimum value for an x-variable. If this number 
is smaller than the number dictated by kinematical limits, the latter is used. The output xlkaleu 
and x2kaleu of the generation routine are the values at which the PDFs should be evaluated. 



Masses and widths 
Kinematics 

Instance of the phase space generator 
Concerns the aeneration of xl,x2 
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B Example of the use of more instances 



As an example of the use of more instances of the program, the user could want to treat several 
subprocesses in one Monte Carlo run, choosing another subprocess for each event. Then the user- 
defined routines could get one more integer input variable iproc determining the subprocess, 
and could look like 



module my_kaleu 
use avh_kaleu 
use avh_kaleu_model 
use avh_kaleu_kinem 
integer , parameter : : r 
type (model_type) , save 
type (kinem_type) , save 
type (kaleu_type) , save 
t ype ( s t r f _t ype ) , s ave 

end module 



c = 20; 
mdl 

kin (nmax) 
obj (nmax) 
str (nmax) 



Masses and widths 
Kinematics 

Instance of the phase space generator 
Concerns the generation of xl,x2 



subroutine my_init ( iproc , process , nfinst , ecm , nbatch, nstep, thrs ) 
use my_kaleu 

integer , intent (in) iproc, process (-2 : 17 ), nfinst , nbatch, nstep 

real (kind ( IdO ) ) , intent (in) :: ecm, thrs 
type (vertx_t ype ) :: vtx 
call my_model ( mdl, vtx ) 

call kaleu_put_process ( mdl, vtx, kin (iproc) , obj (iproc) , process , nfinst , ecm ) 
call kaleu_init_adapt ( mdl, obj (iproc) , nbatch, nstep, thrs ) 
call kaleu_init_strf ( str (iproc) , kin (iproc) , OdO ) 
end subroutine 



subroutine my_gnrt ( iproc , discard , xlkaleu, x2kaleu, pkaleu ) 
use my_kaleu 

integer , intent (in) :: iproc 

logical , intent (out) :: discard 

real (kind ( IdO ) ) , intent (out) :: , xlkaleu, x2kaleu, pkaleu ( : 3 , -2 : 17 ) 
call kaleu_gnrt_strf ( str (iproc) , kin (iproc) , discard , xlkaleu, x2kaleu ) 
call kaleu_gnrt ( mdl, kin (iproc) , obj (iproc) , discard , pkaleu ) 
end subroutine 



subroutine my_wght ( iproc , weight ) 
use my_kaleu 

integer , intent (in) :: iproc 

real (kind ( IdO ) ) , intent (out) :: weight 
real (kind ( IdO ) ) :: ww 

call kaleu_wght ( mdl, kin (iproc) , obj (iproc) , weight ) 
call kaleu_wght_strf ( str (iproc) , ww ) 
weight = weight*ww 
end subroutine 



subroutine my_collect ( iproc , weight ) 
use my_kaleu 

integer , intent (in) :: iproc 

real (kind ( IdO ) ) , intent (in) :: weight 
call kaleu_collect ( obj (iproc) , weight ) 
call kaleu_collect_strf ( str (iproc) , weight ) 
end subroutine 
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